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Abstract 

The form factors parameterizing the B c semileptonic matrix elements can 
be related to a few invariant functions if the decoupling of the spin of the 
heavy quarks in B c and in the mesons produced in the semileptonic decays 
is exploited. We compute the form factors as overlap integral of the meson 
wave-functions obtained using a QCD relativistic potential model, and give 
predictions for semileptonic and non-leptonic B c decay modes. We also discuss 
possible experimental tests of the heavy quark spin symmetry in B c decays. 
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I. INTRODUCTION 



The discovery of the meson by the CDF Collaboration at the Fermilab Tevatron |I| 
opens for interesting investigations concerning the structure of strong and weak interactions 
in the quarkonium-like be hadronic system. The studies will be further developed at the 
hadronic machines currently under construction, such as the LHC accelerator at CERN, 
where a copious production of B c meson and of its radial and orbital excitations is expected 
PJ3]; at these experimental facilities, together with the measurement of the mass of the 
particles belonging to the be (be) family, it will be possible to observe the decay chains 
reaching the 1 So ground state, the B c , which decays weakly. 

A peculiarity of the B c decays, with respect to the decays of the B u ^ and B s mesons, 
is that both the quarks are involved in the weak decay process with analogous probability. 
The weak decays of the charm quark, whose mass is lighter than the b quark mass, are 
mainly governed by the CKM matrix element V ca which is larger than V c b mainly controlling 
the b quark transitions; the result is that both the quark decay processes contribute on a 
comparable footing to the B c decay width. Another peculiar aspect is that the be annihi- 
lation amplitude, proportional to V c b, is enhanced with respect to the analogous amplitude 
describing the B + annihilation mode. 

The above considerations have inspired several theoretical analyses aimed at pre- 
dicting the B c lifetime. Namely, a QCD analysis [0], based on the OPE expansion in the 
inverse mass of the heavy quarks and on the assumption of quark-hadron duality, provides 
for tb c a prediction in agreement (at least within the current experimental accuracy) with 
the CDF measurement: r(B c ) = 0.461^ (stat) ± 0.03 (syst)10~ 12 s 0. The agreement 
supports the overall picture of the inclusive B c decays. 

The calculation of the B c exclusive decay modes can be carried out either using QCD- 
based methods, such as lattice QCD or QCD sum rules, or adopting some constituent quark 
model. So far, lattice QCD has only been employed to calculate the B c purely leptonic width 
||. As for QCD sum rules ||10|| , the B c leptonic constant, as well as the matrix elements 



relevant for the semileptonic decays, were computed in refs. [11,6,12]. These analyses iden- 



tified a difficulty in correctly considering the Coulomb pole contribution in the three-point 
functions needed for the calculation of the semileptonic matrix elements. Attempts aimed at 
taking this correction into account are described in [Q; however, the problem of including 
the contribution of the Coulomb pole for all the values of the squared momentum transfer 
t to the lepton pair has not been solved, yet. Extending to all values of t the expression 
of the Coulomb contribution valid at t max only allows to conclude that it represents a large 
correction to the lowest order quark spectral functions. 

It is worth looking at the outcome of constituent quark models which, although less 
established on the QCD theoretical ground, can nevertheless provide us with significant 
information to be compared to the experimental results. 

The models in refs. [|14l , |T5l have been used in the past to estimate the semileptonic 

B c decay rates. More recently, different versions of the constituent quark model have been 
used to analyze the decays induced both by the b — > c{u) and c — > s(d) transitions |T7| , |T8 



It is noticeable that the calculations can be put on a firmer theoretical ground if some 
dynamical features of the B c decays are taken into account. Such features are mainly 
related to the decoupling of the spin of the heavy quarks of the B c meson, as well as of the 
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meson produced in the semileptonic decays, i.e. mesons belonging to the cc family (r) c , J/i/j, 
etc.) and mesons containing a single heavy quark (B^*\ B^ , £)W). The decoupling occurs 
in the heavy quark limit (m^, m c 3> Aqcd), and produces a symmetry, the heavy quark 
spin symmetry, allowing to relate the form factors governing the B c decays into a 0" and 



1 final meson to a few invariant functions [19]. The main consequence is that the number 



of form factors parameterizing the matrix elements is reduced, and the description of the 
semileptonic transitions is greatly simplified. 

However, at odds of the heavy quark flavour symmetry, holding for heavy-light mesons, 
spin symmetry does not fix the normalization of the form factors at any point of the phase 
space. The normalization, as well as the functional dependence near the zero- recoil point, 
must be computed by some nonperturbative approach. 

So far, the "universal" form factors of semileptonic B c decays have been estimated using 
nonrelativistic meson wave-functions |19l and employing the ISGW model at the zero-recoil 



point [^0|. An analysis in the framework of a different quark model is described in [17]. 

In this paper we present a calculation based on a constituent quark model which has 
been used to describe several aspects of the heavy meson phenomenology The peculiar 
features of the model are related to the interquark potential, which follows general QCD 
properties, such as scalar flavour-independent confinement at large distances, and asymptot- 
ically free QCD coulombic behaviour at short distances. Moreover, the use of the relativistic 
form of the quark kinematics allows to describe heavy-light as well as heavy-heavy mesons, 
and to account for deviations from the nonrelativistic limit. As a result, the B c form factors 
can be written as overlap integrals of meson wave-functions, obtained by solving the wave 
equation defining the model. As discussed in the following, the representation as overlap in- 
tegral of meson wave-functions allows to predict, in the heavy quark limit, the normalization 
of the invariant functions at the zero-recoil point and to obtain, for example, the suppression 
factor between the form factors of the B c transitions into heavy-light mesons with respect 
to the corresponding functions governing the decays B c — > r\ c tv and B c — > J / ifjlv. 

The calculation of the overlap integrals and of the B c semileptonic form factors is pre- 
sented in Sec. Ill, after having reviewed in Sec. II the consequences of the heavy quark spin 
symmetry in B c decays. In Sec. IV, using the obtained invariant functions, we analyze the 
semileptonic decay modes, and in Sec. V, assuming the factorization ansatz, we estimate 
several non-leptonic B c decay rates. Sec. VI is devoted to the conclusions. 



II. HEAVY QUARK SPIN SYMMETRY 

Heavy quark spin symmetry amounts to assume the decoupling between the spin of the 
heavy quarks in the B c meson, since the bc spin-spin interaction vanishes in the infinite 
heavy quark mass limit, as well as the vanishing of the heavy quark-gluon vertex. This 
symmetry has been invoked in Jl9| to work out relations among the semileptonic matrix 
elements between B c and other heavy mesons (both heavy- heavy and heavy-light). The 
main difference with respect to the most well known case of the heavy-light systems is that 
in the latter case one can exploit heavy quark flavour symmetry, which also holds in the 
heavy quark limit and allows to relate B to D form factors. 

In order to apply spin symmetry to B c decays one should distinguish decays due to 
charm transitions from b quark transitions. To the first category belong processes such as 
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B c — > (B s , B*)£v and B c — > (Bd, B^iu, induced at the quark level by the transitions c — > s 
and d, respectively. Since m c <C m^, the energy released in such decays to the final hadronic 
system is much less than m^, and therefore the b quark remains almost unaffected. As a 
consequence, the final B a meson (a is a light SU(3)p index) keeps the same B c four- velocity 
v , apart from a small residual momentum q. The initial and final meson momenta can then 
be written as: p Bc = M Bc v and p Ba = M Ba v + q, with v ■ q = O(-). The relation between 
the residual momentum q and the momentum k transferred to the lepton pair is 

k» = P Bc - P Ba = (M Bc - M Ba K " <f • (2-1) 

In this kinematic situation, exploiting the decoupling of the spin of the heavy quarks in the 
mesons, several relations can be worked out among the semileptonic B c form factors. A 
straightforward way to derive such relations is to use the trace formalism [^,^] Q This has 



been done in ref. [19|, and we repeat here the derivation for the sake of completeness. 

One introduces a 4 x 4 matrix H cb describing the doublet (B c , B*) of cb mesons of four- 



velocity v [0| : 



H c b= y_^ [BT ^_ Bclb] y__ L± ? (22) 

where B*^ and B c annihilate a vector B* and a pseudoscalar B c meson of four-velocity v. 
Under spin rotations of the heavy quarks, H cb transforms as H cb —* S c H cb Sy 

On the other hand, for heavy-light B a and B* mesons, the analogous 4x4 matrix 
describing the (B a , B*) spin multiplet reads: 

H a = { —^ 1 [BT lll - B al5 ] ; (2.3) 



all the fields in ( p.2| ) ,( pT3| ) contain a factor ^jM Bca and have therefore dimension 3/2. 

Applying the trace formalism, one gets that the hadronic matrix elements relative to 
the decays B c — > B)*>lv have the following general form, compatible with heavy quark spin 
symmetry: 



< B a *\v,q\q a Tc\B c ,v>= -yj M B M B Tr[H a QTH cb } (2.4) 

where Q is the most general Dirac matrix proportional to the four-velocity v and to the 
residual momentum q. The calculation using ( ^.2[ ),( P73] ) shows that the various matrix 
elements reduce to: 



< B a ,v,q\V^\B c ,v > = pM Bc 2M B Jp!i v„ + a tt a 2 q„] , 

<B* a ,v, q\V^\B c , v> =-i J2M Be 2M Bi a tt a 2 e^e^gV , (2-5) 



<B* a ,v, qlA^Bc, v > = J2M B 2M B * e* + a M a 2 e* ■ q v,] 



1 For a discussion of the heavy quark formalism applied to the quarkonium system see ref. [24] 
and references therein. 
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where and represent the weak flavour- changing (c — > s, d) vector and axial current, 
respectively, and e is the B* a polarization vector. Therefore, as shown by eq.( [2.5|) , the six 
form factors parameterizing the B c into B a and B* a matrix elements can be expressed in terms 
of two invariant functions, fi" and f^- The main difference with respect to the spin-flavour 
symmetry, holding in heavy-light mesons, is that the normalization of the form factors is 
not predicted at any point of the kinematic range and, in particular, it is not fixed at the 
non-recoil point q = 0. 

Actually, the form factors f2 2 gi ye r i se t° terms proportional to the lepton mass in the 
calculation of the semileptonic rates. Moreover, do not contribute at zero-recoil. The 
scale parameter a is related to the size of the B c meson, it can be assumed as proportional 
to the B c Bohr radius and represents the typical range of variation of the form factors Jl9| . 



The relations (|2.5| ) are valid near the zero-recoil point, where both B c and the meson 
produced in the decay are nearly at rest. In the case of the transitions B c — ► B^*\ the 
physical phase space is quite narrow (the maximum momentum transfer t to the lepton pair 
is t max ~ 1 GeV 2 ) and therefore one can assume that eqs.(|2.5|) completely determine the 
semileptonic matrix elements (modulo a set of corrections mentioned below). The situation 
is different for processes induced, at the quark level, by the b— quark transitions. Let us 
consider the decays B c — > (D, D*)£u, induced by the b — > u transition. In this case, the 
energy released to the final meson is small only near the zero- recoil point, where q 2 <C m 2 c . 
At such kinematic point one can repeat the considerations for the transition B c — > B s £u, 
obtaining the relations: 



< D, v, q\VJB c , v > = v /2M Bc 2M I5 [S 1 v„ + a £ 2 g M ] , 

< D*, v, q\Vfj,\B c , v > = -i J2M Bc 2M D * a S 2 e^e^gV , (2.6) 



< D*, v, qlA^Bc, v > = pM B 2M D , \L X + a S 2 e* ■ q Vfl ] . 

Far from the non-recoil point, the light recoiling quark keeps a large momentum, and there- 
fore terms of the order of — cannot be neglected in the effective theory leading to (|2.6| ) . 

Finally, we consider B c decays into quarkonium states, such as r] c and J/ip. The spin 
decoupling of both the beauty and charm quark allows now to relate the six form factors to 
a single one: 



< r] c ,v,q\V^\B c ,v > = x /2M Bc 2M Vc A 



<J/i/),v,q\A lx \B c ,v>= s /2M Bc 2M J/i) Ael . (2.7) 



Also in this case eqs. (|2.7p are only valid near the zero-recoil point. Nevertheless, in the 
following we use them, as well as eqs. ( |2.6| ), for all physical values of the momentum transfer 
t, in order to compute semileptonic and non-leptonic B c decay rates. This is admittedly a 
strong assumption, and the related uncertainty must be added to the uncertainties coming 
from finite mass and QCD corrections that in principle relate the invariant functions to the 
physical semileptonic matrix elements |]19|| . However, assuming eqs.( |2.7| ) and ( |2.6| ) in the 
whole kinematic range, a number of predictions can be collected; the experimental results 
will then provide us with indications on the numerical importance of the corrections. 
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III. B c FORM FACTORS FROM A CONSTITUENT QUARK MODEL 



In this section we compute the form factors A, and Si by using a relativistic potential 
model which allows to account for two QCD effects. The first one is confinement, which 
produces a suppression, at large distances, of the meson wave-functions, due to the linearly 
increasing interquark potential. The second effect is represented by the deviation of the 
quark dynamics from the nonrelativistic limit. By taking such two effects into account, we 
are able to compute the form factor A in ( |2.7| ) as an overlap integral of B c and J/ip wave- 
functions. Moreover, we can apply the formalism to the transitions B c - flW,5« and 
D£ J at the non- recoil point, and then extrapolate the result to the whole kinematic region 
spanned by the various semileptonic transitions. 

Let us consider A in (|2.7|). In order to compute it, we consider the costituent quark 



model studied in [|2T|, whose essential features can be easily summarized. First, we write 
down an expression for the meson state, in the B+ rest frame, in terms of quark and 
antiquark creation operators, and of a meson wave-function: 



\B+ >=z^^Jdk^ Bc (k) b\-k,r,a) c\k,s,P)\0> (3.1) 



where a and (3 are colour indices, r and s spin indices. The operator b* creates an anti-6 
quark with momentum —k, while c* creates a charm quark with momentum k. A similar 



expression holds for the r\ c (cc) state, as well as for vector 1~ states, as described in |2T 
In the meson state, as written in fl3~T|), the contribution of other Fock states such as, e.g., 
states containing one or more gluons, is neglected. 

The wave-function ip Bc (k) describes the momentum distribution of the quarks in the 
meson. It is obtained by solving the wave equation 



yp + m 2 b + yjfc + ml - M Bc }^ Bc (k) + J dk' V(k, k') ^ Bc {k') = (3.2) 

stemming from the quark-antiquark Bethe-Salpeter equation, in the approximation of an 
istantaneous interaction represented by the potential V . Eq.fl3.2p partially takes into account 
the relativistic behaviour of the quarks in the kinetic term; m c and m^ represent the mass 
of the constituent charm and beauty quark, and M Bc the mass of the bound state. 

The QCD interaction is described assuming a static interquark potential having the form, 
in the coordinate space P3J : 

V(r)= 8n A[Ar-#^l , (3.3) 



33 -2n 



f 



Ar 



with A a scale parameter, n/ the number of active flavours, and the function f(t) given by 

,/ s 4 f°° , sin(qt) r 1 1 i . 

f(t) = - . (3.4) 

it Jo q L m(l + q z ) q Zj 

The interest for this form of the potential is that it continuously interpolates the linearly 
confining behaviour at large distances with the QCD coulombic behaviour at short distances, 
where the logarithmic reduction of the strong coupling constant, due to the asymptotic 
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freedom property of QCD, is implemented. A further smoothing of the potential at short 
distances is adopted, according to quark-hadron duality arguments |fjl . 

The wave equation ( |3.2| ), together with the form ( |3.3| ) of the potential and 



p) of 

the meson state, completely determines the model, which has been extensively studied to 
describe static as well as dynamic properties of mesons containing heavy quarks f26|-|28 



Notice that the spin interaction effects are neglected since, in the case of heavy mesons, 
the chromomagnetic coupling is of the order of the inverse heavy quark masses. Therefore, 
both the pseudoscalar and the vector mesons, being degenerate in mass, are described by 
the same wave-function. 

An equation for the form factor A(q = 0) in ( |2.7| ) can be obtained expressing the b — > c 
flavour- changing weak currents in terms of quark and antiquark operators; for the vector 
current, the expression is 



6, 



a/3 



(27T)* 



dqdq 



m b m c 



E b (q)E c (q> 



[u b (q, r)b\(q, r, a) + v b (q, r)d b (q, r, a)]Y 
[u c (tf, s)b c (if, s, (3) + v c (q, s)dl(q, s, (3)] 



(3.5) 



(E q (k) = ^Jk 2 + rri^, k = \k\); an analogous expression describes the axial current. Then, 
writing down the matrix elements (|2.7|) and applying canonical anticommutation relations 
we obtain: 



21 



A(«r=o) 



i 



2J2M B 2M Vc Jo 



dk 



u Bc (k)u ric (k) {E b + m b )(E c + m c ) - k 2 
jEWc [(^ + m fe )( J B c + m c )] 1 /2 



(3.6) 



where the reduced wave-functions %(lr) are related to the L = wave-functions ipM ac- 
cording to 



u M (k) 



V2% 



(3.7) 



The covariant normalization is adopted: J °° dk\iiM{k)\ 2 = 2Mm- 

The wave- functions ub c and u Vc can be obtained by solving eq. ( |3.2|) by numerical meth- 
ods, choosing the values of the masses m c and m b of the constituent quarks, together with 
the scale parameter A, in such a way that the charmonium and bottomonium spectra are 
reproduced: m b = 4.89 GeV and m c = 1.452 GeV, with A = 397 MeV @ 



A fit of 

the heavy-light meson masses also fixes the values of the constituent light-quark masses: 

115 MeV m 



m 



38 MeV and m. 



It is worth observing that, for the be system, 
all the input parameters needed in ( |3.2|) are fixed from the analysis of other channels, and 
the predictions do not depend on new external quantities. 

The numerical solution of ( |3.2| ) produces the spectrum of the be bound states; the pre- 



-wave resonance are 



2N 



6.28 



dieted mass and the leptonic constant of the first S- 
GeV (the value we use in our analysis) and fg c = 432 MeV, in agreement with other theo- 
retical determinations based on constituent quark models ||29|| , QCD sum rules {M Bc = 6.35 
GeV §) and lattice QCD (M Bc = 6.388 ± 9 ± 98 ± 15 GeV 0). Within the errors" the B c 
mass agrees with the CDF result: M Bc = 6.40 ± 0.39 (stat) ± 0.13 (syst) GeV 0. 
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The obtained B c wave-function UB c {k) is depicted in fig.l. In the same figure we plot the 
wave-functions of the other mesons involved in B c semileptonic decays: B s and Bd, the cc 
states rj c and J /if) together with the first radial excitation v/ c and if)(2S), and the D meson. 

Let us come back to eq. (|3.6| ) which provides the form factor A. For quark masses larger 
than the typical relative quark- ant iquark momentum k, eq.(|3.6|) becomes: 



, 1 / dx V Bc (x) % c (x) , (3.8) 

/2M B 2M Vc J 

where ^ m(x) is defined as 

^m(x) = J dk ifj M (k) . (3.9) 

Eq. ( |3.8| ) shows that the form factor A, at the zero- recoil point, is simply given by the 
overlap integral of the B c and r] c wave-functions in the coordinate space. This result has 



already been obtained in |]T9[] , as it is typical of the calculation of form factors by quark 
models The interest in eq.( |3.8|) is that no factors appear in the integral other 

than the wave functions; this implies that, in the limit where the B c and rj c wave-functions 
are equal (modulo the normalization condition), the form factor A is 1. Although such an 
overlap is not constrained by symmetry arguments, as in the case of the flavour symmetry in 
heavy- light mesons, from eq.( |3.8| ) it turns out that the deviation from unity of the invariant 
function at the zero-recoil point is due to the actual shapes of the meson wave-functions. In 
our specific CcLS6, clS reported in Table |, the deviation from unity is a 5% effect. 

The calculation of A near the zero-recoil point, for a small momentum q, can be per- 
formed by modifying eq. (|3.8|), as discussed in [O: 



A(q) = 1 / dxe^ 2 * Bc (f) *; c (f) , (3.10) 

J2M Bc 2M Vc J 



and using the relation (valid near the zero-recoil point) y = A ^ c ^ c = ^/l + ep/M^. We 
choose to perform an extrapolation of the result in the whole kinematic region, obtaining 
the form factor depicted in fig. 2. The extrapolation provides a form factor having a nearly 
linear (with a small curvature term) y— dependence in the kinematic range of the decays 
B c — > f] c £u and B c — > J/if)£u. 

The same method and the same formulae can be used to calculate the form factor A' of 
B c — > i]' c and B c —>■ ip(2S); the only new ingredient is the wave-function of the ip(2S) radial 
excitation. Due to the oscillating behaviour of u^s), the function A' is suppressed with 
respect to A; interestingly enough, it has a negligible y— dependence, as one can observe in 
fig.2. 

Before discussing the phenomenology of the decays B c — > r] c (J '/ 'if>)£i/ and B c — > 
ri' c (if>(2S))£i>, let us consider the matrix elements relevant for the transitions B c — > B S (B*). A 
feature of the model we are considering is that both heavy-heavy and heavy-light mesons are 



S 



described by the same formalism. Therefore, eg. (|3.6| ) can be applied to calculate = 0), 
substituting m b with m s and the wave-function u Vc with u Bs . In the limit m s — > and for 
a large value of the b— quark mass, eq. (|3.6|) becomes: 



which differs by a factor A= with respect to the analogous relation for A. This factor is a 
consequence of considering a heavy-light meson in the final state instead of a heavy- heavy 
meson, and produces a suppression of the corresponding form factor. Eq.( [3.1lD suggests 
that, for similar (modulo the normalization condition) B c and B s wave-functions, the form 
factor f2f((f = 0) is close to the value f2f((f = 0) = l/v2. The actual value, reported in 
Table |, differs from this value by a 7% effect. 

The two results A(q = 0) ~ 1 and Q{(<f = 0) ~ l/y/2 are the main predictions of our 
analysis. They would deserve independent checks by different theoretical methods, namely 
by QCD sum rules in the heavy quark limit. 

From eq. (|3.11|) it is also possible to derive a relation, proposed in Hl9| , between the form 



factor fif and the leptonic constant of the B s meson. As a matter of fact, in the framework 
of the constituent quark model, the B s leptonic constant, defined by the matrix element: 
< 0|A M |S s (p) >= if Bs Pn, is given by |HJ: 

t f°° ,, , fu ^ (E b + m b )(E a + m s ) ] i/2 fc 2 

IB, = m / dkku Bs {k) — 1-— — r— — r - 3.12 

2nM Bs Jo 1 E b E s J (E b + m b )(E s + m s ) 

For vanishing m s and large m b , f Bs is simply related to the B s wave-function at the origin: 



IB S 



M, 



■*b.(0) , (3.13) 



a relation analogous to the van Royen-Weisskopf formula for the quarkonium state. Ex- 
panding ty Bs {x) near the origin in ( |3.11| ), we obtain: 

1 If 

n{(q = 0) ~ ^=f Bs VM Bs -== J dx V Bc (x) + corrections . (3.14) 



The numerical comparison of ( |3.14| ) with ( |3.11| ), however, suggests that the next-to- leading 



corrections in ( |3.14 ) are sizeable, and therefore the expansion (truncated at the first term) 
leading to eq. fl3.14| ) appears to be of limited usefulness. 

The value of Q{ at zero-recoil is reported in Table |, and the plot of the form factor, 
extrapolated in the whole kinematic region, is depicted in fig.2; the form factor presents a 
soft ^/-dependence in the narrow kinematic range spanned by the semileptonic B c — ► B s , B* 
transitions. 

The same procedure can be applied to compute Qf and Si, and the results are also 
depicted in fig.2. The only new information is that, keeping finite values of the light quark 
masses, a SU(3)p breaking effect between Qf and fl{ of less than 3% is predicted. 

All the invariant functions can be represented by the three-parameter formula 

F{y) = F(0)(l - p\y - 1) + c (y - l) 2 ) (3.15) 
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in terms of the value at zero-recoil, the slope p 2 and the curvature c; the corresponding 
values are collected in Table U 

A remark concerns the invariant functions an d S 2 . As mentioned in Sect. II, such 
form factors do not contribute at the zero- recoil point, since they appear in the term pro- 
portional to the small momentum q. In our approach, based on considering overlap integrals 
of wave-functions of mesons at rest, we cannot provide an independent calculation of fl^ 
and S 2 , which therefore will be neglected in our analysis. Such an approximation, however, 
could have relevant consequences only in the case of the transitions B c — > D^iu; as already 
underlined, for the decays B c — > and B c — > B^ the contribution from f2 2 is always 
proportional to the momentum q, which remains small in these processes. 

Let us conclude the section comparing our form factors A, fll and S x with the outcome 
of the ISGW model |I5| , which has been widely applied to describe the heavy meson decays. 
In the ISGW approach, the form factors exponentially depend on the squared momentum 
transfer to the lepton pair, and at zero-recoil they are given by products of parameters 
relative to the mesons involved in the decays. We depict in fig. 2 the various invariant 
functions obtained in this approach, observing some agreement with our results in the case 
of A; as for the result based on |T!| deviates considerably from the value 1/ y/2 suggested 
by our model. 



IV. B c SEMILEPTONIC DECAYS 

The form factors and flf, A, A'and S x can be used to predict the semileptonic 
B c decay rates, as well as various decay distributions. Before doing the calculation let us 
stress again that an extrapolation is performed for the relevant matrix elements far from 
the symmetry point (zero-recoil) where the form factors are originally computed. Such 
a procedure would require the calculation of the corrections, which could be sizable far 
from the symmetry point, an analysis beyond the aim of the present work. Considering 
the small range of momentum transfer t involved in c — > (s, d) transitions, it is plausible 
that the extrapolation is quite under control for the decays B c — > B^iu, B d *^£u. As for 
B c — > i] c , J/ip£v, the extrapolation is done on a wider range of momentum transfer to the 
lepton pair. However, also in this case it is interesting to make predictions and to compare 
them with the experimental results. Notice that we only consider massless charged leptons 
in the final state. 

Concerning the parameters needed in the analysis, we use the experimental values of 
the masses of 7] c , J/ip, xJj(2S), D(*\ B^ and B s mesons; for the r]' c we use M v > c = 3.66 
GeV, and for M B * we put: M B * = M Ba + (M B * - M Bd ). For the CKM matrix elements we 
use V c b = 0.039 and V u b = 0.0032; the values of V cs and V c d are fixed to V cs = 0.975 and 
Vd = 0.22. The results for the decay widths are reported in Table || where we also report 
the corresponding branching fractions, obtained assuming for r Bc the CDF central value: 
t Bc = 0.46 ps. 

In order to understand the effect of the t— dependence of the form factors, we also report 
in Table |T| the results obtained assuming t— independent invariant functions, with the values 
fixed at the zero- recoil point. The results provide us with an upper bound for the various 
decay widths. As expected, the momentum transfer dependence is mild in the case of the 
B c — > B^*\ B^ decays, where it only provides an effect of less than 10% in the decay rates. 
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This is mainly due to the narrow t— range spanned in such decay modes. In the case of 
B c — ► Vc an d J l^i there is a sizeable effect due to the t— dependence of the form factors. On 
the contrary, in the case of decays into radial excited states, rj' c and i/j(2S), the t dependence 
is negligible. The t— dependence is important for the Cabibbo suppressed B c decays into D 
and D*. 

From Table [II] we conclude that the semileptonic modes are dominated by two channels, 
B c — > B s iv and B c — > B* s lv, in spite of the small phase space available for both the 
transitions; the two modes nearly represent the 60% of the semileptonic width, a result in 
agreement with the predictions available in the literature. 

As for the h — > c induced semileptonic B c transitions, a peculiar role is played by the 
B c decay into J/ 1 ^, due to the clear signature represented by three charged leptons from 
the same decay vertex, two of them coming from J/ijj. This signature has been exploited to 
identify the B c meson at Tevatron |l] , and will be mainly employed at the future colliders [|34| . 



Our prediction for the width of the decay B c — > J/ipiv is: T(B C — > J/if)£u) ~ 21.6 x 1CT 15 
GeV, with an upper bound of 48 x 10~ 15 GeV obtained using a t— independent form factor 
A. The agreement of this result with other calculations in the literature suggests that the 
finite mass corrections, responsible of subleading form factors in the matrix elements, should 
not be large. Tests on the size of such corrections can be performed by measuring the B c 
decay rates into longitudinally and transversely polarized J/ip: T LjT = T(B C — ► J/ipL,T&v), 
together with the corresponding decay distributions. Using the parameterization in 
the decay widths are given by: 



j ^ dy[A{y)f^l[ry-l] 



2 



12tt 3 



G 2 F VlM 5 J/lb rl+S „ , „ 

r T = I dy [A(y)] "J^-^r + 1 - 2 ( 41 ) 

where r = M Bc /Mj/^ and 5 = WJ^ ■ ^ ne measurement of dTi/dy provides infor- 

mation on A and V c b] in particular, if the curvature term in A(y) is neglected, the ratio 
Tt/Tl gives access to the slope p 2 . The combination V c bA(l) can be obtained from the 
measurement of Tl and from the total width, and therefore a measurement of V c b is possible 
using this decay channel [Q,^]. Such new determinations of the CKM element V c b, even 



though not accurate as from Bd and B u decays, would represent an important consistency 
check of the Standard Model. 

Tests of the spin symmetry are provided by the measurement of the decay distributions 
in the y variable, whose deviations from the distributions related to a unique form factor A 
would imply the presence of spin symmetry-breaking terms. 

Let us finally observe that our prediction for the rates of the decays into 0~ (cc) states, 
B c — > rj c lv and B c — ► T]' c £u, is smaller than the value reported by other analyses. 



V. NON-LEPTONIC B c DECAYS 

Estimates of the decay rates of several two-body non-leptonic B c transitions can be 
obtained adopting the factorization approximation. Such an approximation finds theoretical 
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support in few cases (large N c limit; mj — > oo limit in 6 — > it transitions involving heavy-light 
meson systems H35|); nevertheless, it is widely used to estimate non-leptonic decay rates of 



mesons containing heavy quarks. 

Let us first consider non-leptonic B c decay modes induced, at the quark level, by the 
b —>■ c and u transitions. The effective Hamiltonian governing the processes reads: 

H eff = ^ { VcbiaMQ? + c 2 (fi)Qf] + V uh [ Cl {v)Qf + c 2 {ix)Qf) + h.c.} 

+ penguin operators ; (5-1) 

Gp is the Fermi constant, are CKM matrix elements and Cj(/i) scale-dependent Wilson 
coefficients. The four-quark operators Qf 1 and Q% are given by 

Ql = K d {du)y. A + {SU) V „ A + {dc)y. A + V c * s (Sc) V _ A ] {cb)y_ A 

Qf = Kd (cu) V -A (db) v ^ A + V* s (cu) V -a (sb)v„ A + V c * d (cc) V -a (db) + V c * s (cc) V -a (sb)] 

(5.2) 

with {q\q 2 )v-A = <?i7^(l — 15)12] analogous relations hold for Qf' and Qf 3 . 

As well known, the factorization approximation amounts to evaluate the matrix elements 
of the four-quark operators in ( |5.2j ) between the initial B c state and the final two-body 
hadronic states as the product of quark-current matrix elements. We adopt this approxima- 
tion in the calculation of the rates, neglecting the contribution of penguin operators, since 
their Wilson coefficients are small with respect to c\ and c 2 (interference effects of penguin 
diagrams are of prime importance in producing CP violating asymmetries in B c decays). 
Moreover, we do not take into account the weak annihilation contribution represented by 
a B c meson annihilating into a charged W; in this amplitude, the final hadronic state is 
entirely produced out of the vacuum, and therefore the contribution should be character- 
ized by a sizeable form factor suppression. Annihilation processes are presumably relevant 
mainly for rare or suppressed B c decays; in these cases they deserve a dedicated analysis. 

A further remark concerns the Wilson coefficients ci(/i) and c 2 {^)- Writing the factorized 
amplitudes and taking into account the contribution of the Fierz reordered currents, it turns 
out that the relevant coefficients are the combinations: ai = C\ +C,c 2 and a 2 = c 2 + £ci, with 
the QCD parameter £ given by £ = 1/N C . Several discussions concerning this parameter are 
available in the literature. We choose a>i = C\ and a 2 = c 2 , i.e. £ = 0, in the spirit of the 
large iV c limit, and use c\ and c 2 computed at an energy scale of the order of m^. A detailed 
analysis of 1/N C corrections to the coefficients ai,a 2 as well as of the role of color-octet 



current operators in B decays can be found in [36|. Analogous considerations hold for the 
decays induced by the c — > s(d) transitions; in this case we choose the coefficients c\ and c 2 
at the scale of the charm mass. 

The factorized amplitudes can be expressed in terms of the form factors in eqs.( |2.5| ), ( |2.6| ) 
and ( |2.7|) , and of leptonic decay constants defined by the matrix elements: < 0|v4 M |M(p) >= 
ifuVp and < 0|T^|V(p,e) >= f v M v e^. We use the values: f w + = 0.131 GeV, f p + = 0.208 
GeV and f ai = 0.229 GeV; f K+ = 0.159 GeV, f K *+ = 0.214 GeV and f Kl = 0.229 GeV; 
f Vc = 0.31 GeV, f rfc = 0.23 GeV, = 0.38 GeV, = 0.28 GeV, and finally f D = 0.2 
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GeV, f Ds = 0.24 GeV and f D » = 0.23 GeV, f D * = 0.275 GeV. Such values correspond to 
experimental results or to average values from lattice QCD and QCD sum rules 0. 

The decay rates of several non-leptonic B c transitions, obtained using ci(m&) = 1.132, 
Cii^b) = —0.286 and Ci(m c ) = 1.351, c 2 (m c ) = —0.631, are collected in Tables |T|, [TV- 



Also in this case we use the physical phase space together with the expression of the matrix 
elements in 

Few comments are in order. We observe the dominance of the decay modes induced by the 
charm transition, and in particular of the channel B^ — > B*p + , which represents more than 
10% of the total B c width. It would be interesting to experimentally confirm this prediction, 
even though the final state presents severe reconstruction difficulties. From the experimental 
point of view, more promising are the decay modes having a J/ip meson in the final state; 
among such modes, the decay channels B^ — > J/ip7r + and B+ — > J/ipP + are particularly 
useful for the precise measurement of the B c mass, by the complete reconstruction of the 
final state. Also the decay into a x is of particular interest, due to the large decay rate. 

Several tests of factorization can be carried out, mainly using the decay channels having 
a J/ip in the final state. For example, the assumption of the factorization approximation, 
together with the heavy quark spin symmetry, implies that the relation 

T(B+ -> J/^+) _ ^V 2 d a\fl (5 3) 



dT{BZ^JH>l+v) | M Bc M J/V , 

dy ' y—VT c 

M 2 +M 2 

holds in the limit M n — > (y n = 2 ®i B W^ )- ^n ana l°gous relation holds for the B c 
transition into the radial excited state ip(2S): 

Y{Bt - ^(2S)tt+) ^V* d a\fl 



dT(B+^(2S)l+i/) j M Bc M^(2S) 



(5.4) 



In the case of a p meson in the final state one has: 

T(B+ -> J/^p+) _ 37rV>?/ p 2 [8M^M p 2 + (M 2 c - M 2 J/f - M 



I?} 



dv{B+^j/w+v) | smI M 5 T/ , 

d y \y=y P Uc J lv 



\Jy 2 — l[r 2 y 2 — 6ry p + 2r 2 + 3] 



(5.5) 



A being the triangular function, r = 1 f, Bc and y a = ^."t 3 Jf — -. 

To test eqs.(|5~3| ) -( |5.5|) two-body decay rates and the differential B+ — > J/i[)£ + v decay 
width are required; the measurement of such quantities, possible at the hadronic facilities, 
would provide us with important information on the heavy quark spin symmetry as well as 
on the factorization approximation in B c decays. 



2 A description of the current theoretical situation concerning the heavy meson leptonic decay 



constants is reported in the Appendices C and D of ref. [33]. 
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VI. CONCLUSIONS 



We have presented a determination of the invariant functions parameterizing the semilep- 
tonic B c matrix elements in the infinite heavy quark mass limit. The form factors are ob- 
tained as overlap integrals of meson wave-functions, obtained in the framework of a QCD rel- 
ativistic potential model. An interesting result is that, although not constrained by symme- 
try arguments, the normalization of the form factor A describing the transition B c — > J / ipii/ 
is close to 1 at the zero- recoil point, as being the overlap of similar wave-functions. On the 
contrary, the form factors relative to the transitions into heavy-light mesons, at zero-recoil 
point, are suppressed by a factor ~ l/\/2 with respect to A. These results have several phe- 
nomenological consequences, in semileptonic and non-leptonic B c decay processes, which can 
be experimentally tested. Moreover, they affect other important processes, such as radiative 
flavour- changing B c decays |37|] and CP violating B c transitions [|38| , jl8|l . In particular, the 
invariant functions computed in this paper can be useful to identify the B c decay channels 
characterized by a clean experimental signature, a large branching fraction and a visible CP 
asymmetry; the identification of this kind of decay modes is of paramount importance for 
the physics program of the experiments at the future accelerators. 
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TABLES 



TABLE I. Parameters of the form factors (ip' = ip(2S)). The functional dependence is in ( |3.15| ). 



Channel 


Form factor 


F(l) 


P 2 


c 


B c -» B S (B* S ) 


Of 


0.66 


8 





B c -> B d {B* d ) 


nf 


0.66 


8 





B c -» Vc(J/^) 


A 


0.94 


2.9 


3 


Be - 7/ c (V) 


A' 


0.23 








£ c -> L>(L»*) 


Si 


0.59 


1.3 


0.4 



TABLE II. Semileptonic Bj decay widths and branching fractions. 



Channel 


L(10- 15 GeV) 


L L (10- 15 GeV) 


L T (10- 15 GeV) 


BR 


B+ -> B s e + v 


11.1(12.9) 






0.8(0.9) x 10" 2 


B + - fije+x, 


33.5(37.0) 


19.1(21.4) 


7.2(7.8) 


2.3(2.5) x 10~ 2 


£+ -» 5 d e+i/ 


0.9(1.0) 






0.06(0.07) x 10~ 2 




2.8(3.2) 


1.6(1.8) 


0.6(0.8) 


0.19(0.22) x 10~ 2 


B+ -> T/ c e+z^ 


2.1(6.9) 






0.15(0.5) x 10" 2 


-> J/^e+^ 


21.6(48.3) 


13.2(33.2) 


4.2(7.6) 


1.5(3.3) x 10~ 2 


S+ -> r£e+i/ 


0.3(0.3) 






0.02(0.02) x 10~ 2 


5+ -> ^'e+i/ 


1.7(1.7) 


1.1(1.1) 


0.3(0.3) 


0.12(0.12) x 10~ 2 


5+ -> D°e+^ 


0.005(0.03) 






0.0003(0.002) x 10~ 2 


5+ -> £>*°e + z/ 


0.12(0.5) 


0.08(0.35) 


0.02(0.05) 


0.008(0.03) x 10~ 2 
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TABLE III. Non-leptonic (b — ► c, u) decay widths and branching fractions. 



Channel 


r(10~ 15 GeV) 


BR 


Channel 


r(10" 15 GeV) 


BR 




i] c ir + 

VcP + 

Vcaf 


, '2 n 90 

af 0.75 
af 0.96 


2.6 x 10" 

6.7 x 10" 
8.6 x 10" 


-4 
-4 
-4 




(2^ U.UZo 

af 0.041 
af 0.05 


O V 1 n~ 
Z X 1U 

3.6 x 10 
4.4 x 10 


5 

-5 
-5 


'Ic 71 

VcP + 

Vc a l 


n 2 D 074 
af 0.16 
af 0.15 


6.6 x 10" 
1.5 x 10" 
1.4 x 10" 


-5 
-4 
-4 


n' K + 

r/ c K*+ 


a 2 OD^ 

til U.UUiJU 

af 0.008 
af 0.0075 


7.4 x 10 
6.7 x 10 


6 

-6 

-6 


T /„/,_ + 

J/ipaf 


n 2 1 zLS 

af 4.14 
af 5.78 


1.3 x 10" 
3.7 x 10" 
5.2 x 10" 


-3 
-3 
-3 


77*/] ?T+ 

j/i>K*+ 


n 2 n n7fi 
a-^ u.uiu 

af 0.23 

af 0.3 


O.O X 1U 

2 x 10" 
2.7 x 10 


-5 
4 

-4 


W P 
Vat 


af 0.22 
a^ u.04 
af 0.65 


1.9 x 10" 
4.86 x 10 
5.8 x 10" 


-4 
-4 
-4 


ii'Kt 


af 0.01 

n 2 n (YX 

a^ u.uo 
af 0.033 


9.3 x 10 
9 ft v in 

3 x 10" 


-6 
-5 
5 


D+D° 
D + D*° 
r>*+ h° 


a 2 . 0.15 
a 2 , 0.13 

n 2 1 4fi 

al 2.4 


8.4 x 10" 

7.5 x 10" 
8.4 x 10" 
51.4 x 10 


-6 
-6 
-5 
-4 


DfD° 
D+D*° 
r>*+ f)° 

D*+D*° 


a 2 2 0.01 
a\ 0.009 
n 2 n ns7 

u-2 U.UO / 

al 0.15 


6 x 10" 

5.3 x 10 

t; v in- 

O A ±U 

8.4 x 10 


7 

-7 
6 

-6 


Vc D s 


(ai 7.8 + a 2 1.6) 2 x 10" 1 
(ai 3.6 + a 2 6.05) 2 x 10" 1 
(ai 1.5 + a 2 o.zj x 10 
(ai 0.79 + a 2 1.8) 2 x 10" 1 


5 x 10 _ 
3.8 x 10" 
3.7 x 10" 

1 x 10- 


3 

-4 
-5 

5 


Vc D+ 

Vc D 

V ' C D*+ 


(ai 0.86 + a 2 0.46) 2 x 10" 1 
(ai 0.7 + a 2 0.9) 2 x 10" 1 
(ai 0.28 + a 2 0.7) x 10 
(ai 0.17 + a 2 0.8) 2 x 10" 1 


5 x 10- 
2 x 10" 

1 X 10 

6 x 10- 


5 
5 
6 
8 


J/^D S 
J/^jD* 
$D S 


(ai 6.7 + a 2 2.3) 2 x 10" 1 
(ai 11 + a 2 10.4) 2 x 10" 1 
(ai 1.4 + a 2 1.33) 2 x 10" 1 
(ai 2.75 + a 2 7.8) 2 x 10" 1 


3.4 x 10" 
5.9 x 10" 
1 x 10- 
5.7 x 10" 


-3 
-3 
4 

-5 


J/4>D+ 
J/ipD*+ 

t(j'D*+ 


(ai 1.31 + 02 0.47) 2 x 10" 1 
(ai 2.02 + a 2 2.3) 2 x 10" 1 
(ai 0.35 + a 2 0.36) 2 x 10" 1 
(ai 0.55 + a 2 1.76) 2 x 10" 1 


1.3 x 10 
1.9 x 10 
5.8 x 10 
8.7 x 10 


-4 
-4 
-6 

-7 



TABLE IV. Non-leptonic (c — ► s,d) B+ decay widths and branching fractions. 



Channel 


r(10~ 15 GeV) 


BR 


Channel 


r(10~ 15 GeV) 


BR 


B s ir + 


a\ 30.6 


4 x 10~ 2 


B S K+ 


a\ 2.15 


2.7 x 10~ 3 


B s p + 


a\ 13.6 


1.7 x 10~ 2 


B S K*+ 


a\ 0.043 


5.4 x 10~ 5 




a\ 35.6 


4.5 x 10 


B*K + 


a\ 1.6 


2 x 10~ 3 


B* sP + 


a\ 110.1 


1.4 x 10" 1 








B d 7T+ 


a? 1.97 


2.5 x 10" 3 


B d K+ 


a\ 0.14 


1.8 x 10" 4 


B dP + 


a\ 1.54 


2 x 10~ 3 


B d K* + 


a? 0.032 


4 x 10~ 5 


Bp + 


af 2.4 


3 x 10~ 3 


B d K + 


a 2 0.12 


1.6 x 10" 4 


B* d P + 


a\ 8.6 


1 x lO- 2 


B*K*+ 


a\ 0.34 


4.4 x 10- 4 
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FIGURE CAPTIONS 



Fig. 1 

Reduced L = wave-functions %(fc) of heavy-heavy (B c , J/^, ip(2S)) and heavy-light (B s , 
Bd, D) mesons. The wave-functions are obtained by solving the wave equation ( p.2|) ; they 
describe both the pseudoscalar 0~ and vector 1~ mesons. 



Fig. 2 

Form factors of B c semileptonic decays. The variable y is related to the squared momentum 



t, transferred to the lepton pair, by the relation: y 



M l+ M lr 



The solid lines correspond 



2M Bc M M 

to the form factors obtained by the model discussed in the paper; the dashed lines refer to 
the model in ref. |15 
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